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Abstract 

The Hubbard model has been employed successfully to understand many aspects of 
correlation driven physical properties, in particular, the magnetic order in itenerant elec- 
tron systems. In some systems such as Heusler alloys, manganites etc., it is known that, 
in addition to magnetic order, distortion induced by Jahn- Teller (J-T) effect also exists. 
In this paper, based on two-fold degenerate Hubbard model, the influence of magnetic 
order on J-T distortion is investigated. The electron correlation is treated using a spec- 
tral density approach and J-T interaction is added to the model. We find that magnetic 
order and structural distortion coexist at low temperature T for a certain range of electron 
correlation strength U, J-T coupling strength G and band occupation n. At T = 0, for 
a given n and U, magnetic order is present but distortion appears only for a G larger 
than a critical value. We also studied the temperature dependence of lattice strain and 
magnetization choosing a G close to the critical value. 
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1 Introduction 



Strongly correlated electron systems have been at the centre of condensed matter research for 
quite a long time. One of the oldest models to describe correlations among itenerant electrons 
is the Hubbard model. This model is very well researched using practically all the available 
theoretical techniques of condensed matter. It is proved to be relevant in many contexts such 
as itenerant magnetism, metal insulator transitions, high temperature superconductivity etc. 
On the other hand, the interplay of magnetic and electric phenomena has attracted attention 
in the context of colossal magneto resistance (CMR), the prototype materials for this being the 
manganites[l,2]. In these materials it is still not quite clear whether the observed magnetism 
stems from localized moments (t2 9 -electrons) or from the itenerant ones (e s -electrons). In 
manganites another essential interaction is the Jahn- Teller (J-T) coupling[3] which causes a 
structural distortion. The story of manganites is therefore, not complete without including 
magnetism and J-T coupling. If the magnetism is due to the band electrons and if the band 
degeneracy can be lifted by J-T coupling, then the study of Hubbard model extended by a J-T 
term is certainly relevant to these systems. This model is interesting from another point of 
view also. The emerging field of orbitronics, like CMR, also has potential applications. The 
orbital dynamics of active orbitals is strongly influenced by their interaction with lattice. One 
way of including electron-lattice interaction is via J-T coupling [4-7]. Here again to be realistic, 
one has to take the correlation among the orbitals into account, which brings us back to the 
Hubbard model extended by J-T interaction. Heusler alloys [8, 9] are another class of systems 
where there is a strong interplay between magnetism and lattice distortion. Here of course, the 
magnetism is that of a local moment system. In this context also, the general effect of magnetic 
order on lattice distortion is worth studying. 

The Hubbard model in general is obviously not exactly solvable. One of the approximations 
proposed is the spectral density approach (SDA) [10]. In this approach, an ansatz is made 
that the spectral density is predominantly a two-peak function which is a very reliable ansatz 
[10,11]. The limitation of this ansatz is that, in this approximation, the quasiparticles are 
infinitely long living. Earlier investigations have shown that the lifetime of quasiparticles does 
not qualitatively influence magnetic ordering [12,13]. The J-T coupling can be incorporated 
into the Hubbard model by renormalizing the band energies. In the present work, the SDA 
is used to solve the renormalized Hubbard model. The influence of the correlation driven 
magnetic order on the structural distortion caused by J-T interaction is studied by calculating 
the magnetization and strain self-consistently. 



2 The model and theory 



The model consists of correlated electrons in a doubly degenerate band interacting with the 
lattice via the J-T interaction, 

H = H S + H JT + H L (1) 
where H s is the Hubbard Hamiltonian for two degenerate bands and is given by 

a,i,j,a a,i,a 

is the hopping integral for hopping of the electrons from lattice site % to j. c^^Caia) is the 
creation(annihilation) operator for an electron in the a-state on the lattice site i with spin a. 
a — 1, 2 is the band index, /i is the chemical potential. n a i a = ^ a i a c a i a is the number operator. 
U is the intra-atomic Coulomb repulsion which we take to be the same for both the bands. The 
intersite hopping integrals are connected to Bloch energies e(k) by 

e(k) = ^E r «e-*- (Ri - Rj) (3) 

i,3 

The electron density in the degenerate band couples to the static elastic strain through the J-T 
interaction. In the case of a tetragonal distortion, this interaction is described by [5- 7] 

H JT = Ge^2(n lia - n 2i(T ) = Ge^(ni k(J - n 2k a)- (4) 

i,cr k,(T 

G is the strength of the J-T coupling and e is the lattice strain given by 

G 

e = ^77T E(< n ^ > ~ < n 2i<x >) (5) 

where C is the elastic constant which we take to be unity for numerical calculations. It is 
clear that H JT tries to create a difference in the occupation of the two degenerate bands. The 
difference in occupation leads to the building up of the strain. Thus, under suitable conditions, 
there is a spontaneous splitting of the bands and building up of the strain which indicates a 
structural transition. The building up of the strain, however, leads to an increase in the lattice 
elastic energy [4-7] which is given by 

H L = ±NC e 2 (6) 

Where N is the total number of atoms. Since this term is a c-number and we are not looking for 
the ground state whose energy has to be minimum, we leave this term out of our consideration. 
Without resorting to any approximation, we can absorb the J-T term into H s by modifying the 
band energies for the two bands as 

e a (k) = e(k) + {-l) a Ge. (7) 



Where Ge, is half the separation between the centres of the gravity of the split free-bands. 
With this, the model Hamiltonian becomes, 

H = E ( e «( k ) - /")cika C ak CT + ^ U Y1 n ai*n ai „ a , (8) 

which is a two-band Hubbard model. The two bands are coupled through the strain which 
is introduced J-T by interaction. In order to calculate the strain one has to calculate the 
single-electron Greens function 

G a ka(E) = ((c aka ; cl ka )) E = — — ■ (9) 

h, - e a {k) + fx — 2^aa{ti) 

That means one has to calculate the self-energy T, aa (E). In the recent past several approxima- 
tion schemes have been developed to calculate self-energy. We make use of the self-energy which 
is obtained[5] by means of the spectral density approach(SDA). As far as magnetic properties 
are concerned, this approach has proven to be quite reliable. 

Z aa (E) = Un a -„ + ^V-.(l-na- g ) 

E + fi- B a _ a -U{1- n a _ a ) 

The pioneering work of Harris and Lange [11] to the Hubbard model demonstrates that, in 
the strong coupling regime, the single-electron spectral density consists of two prominent peaks 
near the energies T and T + U . This is the starting point for the SDA which uses a linear 
combination of two weighted ^-functions. Use of 5-functions implies neglect of quasiparticle 
damping. This is not a serious restriction, at least for strongly coupled systems. It is known 
that the magnetic stability is adversely affected by finite quasiparticle lifetimes quantitatively 
but the qualitative behaviour remains the same. The ansatz contains four unknown parameters, 
namely, the two quasiparticle energies and the two spectral weights, which are fitted by equating 
the first four, rigorously calculated spectral moments. It turns out that the self-energy given 
above (Eq.(lO)) fulfills the correct high-energy behaviour [12]. The SDA is essentially equivalent 
to the Roth two-pole approximation for the Greens function [14,15]. This SDA self-energy can 
also be obtained by using the Mori-Zwanzig projection technique [16-18]. B a - a is known as the 
so-called spin-dependent band shift which is decisive for ferromagnetism and is given by 

1 f+OO 



Baa T aa 



im aa (l - n aa ) 



3 / °°dEf_(E)\-E aa (E)-l 
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(E-T aa + pi- E aa (E))G a(7 (E) - l] (11) 



Where T aa = T + (— l) a Ge. T is the centre of gravity of free Bloch band. G aa (E) = 
jr J2kG a kcr(E). From G a - k(T (E) one can obtain the density of states by using 

Paff (E) = -4f E 3G akCT (£ + t0+) (12) 

7TiV , 



From the knowledge of the density of states, the expectation values can be evaluated: 

(rw) = J dEf_(E) Paa (E) (13) 

Where f-(E) = 1/ (l + e^ E ^j is the Fermi function with f3 = 1/kT. The chemical potential /i 
is fixed by the constraint 

n = E (rw) = constant. (14) 
Now, in terms of the average occupations of the bands (rw), the strain e is given by 

e = -^E«^>-<^» (15) 

The magnetization m is the sum of the magnetizations m x and m 2 of the split bands which are 
defined by 

m a = (n at ) - (n ai ) (16) 

For the numerical evaluation of orbital occupancies and band shift, the k-summation can be 
conveniently replaced by an integration over energy using the model density of states: 

Po(£) = ^£<H£-e(k)) (17) 
We solve Eqns.(14), (15) and (11) self-consistently. In the next section we present the results. 



3 Results and discussion 

The aim of the present investigation is to study the influence of magnetic order on the structural 
transition. The model system is characterized by the following parameters: the density of states 
for free Bloch band which is chosen to be of bcc type [19] with a width W — 1, the intra atomic 
Coulomb repulsion U, electron-phonon coupling constant G and the Band filling n. For various 
sets of these model parameters we have evaluated the band occupancies and the strain self- 
consistently. 

First we consider the case of G = i.e the case of no J-T splitting to find the parameter 
space in which the system is magnetically ordered. In the absence of J-T distortion, the bands 
are degenerate and both the bands have same occupations. We determined the T c which is the 
temperature at which the spontaneous magnetization becomes zero. In Fig.l we have plotted 
Tc as a function of intra-atomic Coulomb repulsion U for various values of band filling n. It 
is noticed that for a given n it requires a minimum U for the onset of magnetic order and the 
T c increases with increase of U and gets saturated for a large value of U . From the inset of 
Fig.l it can be noted that only a certain range of the values of charge carrier concentrations 
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Figure 1: Curie Temperature Tc as a function of Coulomb interaction U for different values 
of band filling n without Jahn- Teller distortion(i.e. with G — 0). The inset shows the n 
dependence of T c for U = 3. 

is favourable for magnetism. This can be understood from the T c dependence of n for U = 3 
which is plotted in the inset. 

In order to see the influence of magnetic order on the structural transition we start with the 
situation where the model system is magnetically ordered and then switch on the the electron- 
phonon coupling constant G which makes the system to be eventually distorted. For this we 
choose the model system with n = 1.4 and U = 3 where the T c = 290 for G = and calculate 
e self-consistently from Eq.(15) by progressively increasing the value of G at T = 0. The result 
is depicted in Fig. 2 (upper part). It is found that J-T splitting does not take place unless the 
electron-phonon coupling is larger than a critical value G c (~ 0.462). When the value of G is 
further increased, the strain increases almost linearly. Below the critical value G c both the 
bands are magnetically ordered with equal magnetization(mi = vri2 ^ and e = 0). However, 
above the critical value of G c the magnetization is mostly due to the first band and is negligible 
for the second band(mi 7^ 0, m 2 ps and e 7^ 0). With further increase in G the magnetization 
of the first band increases slowly while the magnetization of second band approaches zero, and 
the strain exhibits near linear growth. These are all due to the redistribution of electrons from 
the second band to the first band. All these features can be understood from the qusiparticle 
density of states plotted in Fig. 3. The topmost part of DOS which corresponds to G < G c shows 
that the DOS is the same for both bands(l and 2) for a given spin direction. The spectral weights 
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Figure 2: Strain and magnetization as functions of electron-phonon coupling constant G for 
band filling n = 1.4 and U = 3 at T = 0. 

are, however, different for the two spin directions so that we get a non-zero magnetization for 
both the bands and at the same time, the strain is zero due to equal population of bands 1 
and 2. The position of the upper Hubbard bands is determined by magnitude of U and since 
U is large in our case, they are never occupied at T = but the spectral weights of the two 
bands are determined by the sum rule. In this sense the upper Hubbard bands influence the 
physics of the problem. For example the DOS of spin -up in upper Hubbard band is nearly 
a 5-function whose weight depends on the electron concentration. For G > G c the band 1 is 
more populated than the band 2 leading to a finite strain. The spin-up and spin-down parts of 
the band 1 are unequally filled leading to a finite magnetization m\. However, for the band 2, 
the DOS for spin-up and spin-down are the same (dotted line for the spin-down overlaps with 
full-line for the spin- up case) and this makes rri2 = 0. 

From the study of the G dependence of the strain at T = (Fig.2), it can be concluded 
that the presence of magnetization causes a redistribution of electrons between the bands by 
creating a population difference between the spin levels. This redistribution suppresses the 
strain. In order to study the interplay between structural and magnetic transitions, we study 
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Figure 3: Quasiparticle density of states(for the lower subband in the positive half of the frame 
and for the upper subband in the negative half of the frame) as a function of energy for various 
values of electron-phonon coupling constant G. Full lines for spin up and dotted line for spin 
down bands. Thin line is the Fermi function, n = 1.4, U = 3 and T = 



the temperature dependence of the strain and the spontaneous magnetization for different 
values of G, namely (i) G > G c , (ii) G > G c but close to G c and (iii) G < G c but close to G c 
for U = 3 and n = 1.4. We note that G c = 0.462 for the above parameters. For G > G c , at 
T = 0, as pointed out earlier, m! ^ 0,m 2 ^ and e ^ 0. With the increase of temperature, 
both mi and e decrease (Fig. 4a) slowly at first because of redistribution electrons among spin 
and orbital levels of bands. At T c , the magnetization falls to zero abruptly and e rises in a 
discontinuous manner. It is interesting to note that in the absence of J-T interaction, the 
magnetization vanishes smoothly at a higher value of T c (dotted line in right panel of Fig.4a). 
The J-T interaction suppresses T c and also the magnitude of magnetization which indicates 
that the two order parameters compete. Above T c , the strain is high and vanishes smoothly 
at T s , the transition temperature from distorted to undistorted phases. For < T < T c there 
is coexistence of magnetic order and lattice distortion. As T T c the free energy of the state 
with lattice distortion is lowered than that of coexistence phase. We note that the coexistence 
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Figure 4: Temperature dependence of strain and magnetization for various values of G. n = 1.4 
and U = 3. a) The dotted line is the magnetization in the absence of J-T distortion. G = 
0.49. b) The dotted and dashed lines are the magnetization of sub bands land 2 respectively. 
G = 0.47. c) Both the sub bands of equal magnetization which is denoted by the dotted line. 
G = 0.45. 



of magnetic order and lattice distortion has been observed in La\- x Ca x MnO^ at x — 0.15 
[20]. Similar results have been reported [21] recently for the two-band Kondo-lattice model by 
using an interpolating self-energy approach with a modified Ruderman-Kittel-Kasuya-Yosida 
mechanism. 

When G > G c but close to G c , the behaviour of e and m are different (Fig.4b) compared 
to for G > G c at finite but low temperature. To start with, the magnetization mi ^ 0, but 
m2 ~ and e ^ 0. However, e is smaller than in the case of G > G c at low T. When the 
temperature reaches a certain value, e abruptly drops down, mi decreases slightly and m 2 
becomes substantially different from zero so that the total magnetization increases sharply. 
With further increase of T, the magnetization disappears at T c in a discontinuous fashion 
and e is also augmented in a similar way. Thus, it is found that a thermal enhancement of 
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Figure 5: Quasiparticle density of states(for the lower subband in the positive half of the frame 
and for the upper subband in the negative half of the frame) as a function of energy at various 
temperatures for G = 0.47. Full lines for spin up and dotted lines for spin down. Thin line is 
the Fermi function, n = 1.4, U — 3 and T = 0. 



magnetization can exist within a certain range of temperature. This interplay of magnetization 
and distortion is intimately related to the non-linear dependence of spin -level splitting on 
the J-T-splitting. In other words, the J-T interaction renormalizes the correlation in such a 
way that the correlation decreases in band 2 which is shifted to higher energy due to the J-T 
interaction. In such a situation, the magnetization may disappear in this band. However, when 
the J-T-splitting becomes smaller at higher T, the effective correlation in the band 2 goes up 
leading to a finite magnetization. These arguments are supported by the thermal behaviour 
of the DOS (Fig.5). At low T (T = 50) the states with both the spin directions of band 2 
are equally populated and therefore m<i = 0. On the other hand, in case of the band 1, there 
is unequal population of the spin sublevels as is evident from the position of the Fermi edge. 
As T increases, for example at T = 150, there is a relative shift of the spin-down DOS of the 
band 2 so that m-i become finite. There is also a small decrease in m\. Again at T = 250, the 
DOS for both the spin directions is the same so that the total magnetization m = since the 



spin-down population has increased compared to what it was at T = 50. 

For the case G < G c but close to G c (Fig. 4c), both the bands are magnetic and the strain 
is zero for < T < T c . Above T c , the strain appears in a discontinuous fashion. When G 
is decreased further, the strain is always zero for the whole temperature region, even though 
G ^ 0. The re-entrant behaviour of e shown in Fig.4c is observable only in a narrow range 
of G values close to G c . When G is below this range, the magnetic state is uneffected by J-T 
interaction. This demonstrate that the coexistence and mutual competition exist only in a 
narrow range of parameter space. 

4 Conclusions 

The J-T interaction is added to the Hubbard model to study the effect of magnetic order on 
structural distortion in a strongly correlated electron system. The J-T interaction is absorbed 
into the Hubbard model by renormalizing the band energies. Then the Hubbard model is solved 
using the spectral density approximation. Using this approximation, the spontaneous strain and 
magnetization are calculated selfconsistently for various sets of model parameters. The results 
reported earlier, namely, that ferromagnetic order requires a critical Coulomb repulsion U and 
is present only for a certain range of band filling [10], and that there is a spontaneous splitting 
of the degenerate bands only for a J-T coupling G larger than a critical value [21] are recovered. 
In the strong coupling limit and for small values of G the magnetic state is unaltered in a narrow 
range of G for a given value of U and n and there is a competition between magnetic order and 
J-T distortion when they coexist. This is manifested in a reduction of T c and magnetization. 
Within a finite temperature interval an enhancement of magnetization in the magnetic state is 
found. The general features of these results are in tune with the experimental results in doped 
manganites [20]. 
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